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We compute the O (Qs) corrections to the differential rate of the semileptonic decay 6 clvg. at 
the "intermediate recoil" point, where the c-quark mass and the invariant mass of the leptons are 
equal. The calculation is based on an expansion around two opposite limits of the quark masses 
m-6,c: rric — mt and nic <^ mi,. The former case was previously studied; we correct and extend that 
result. The latter case is new. The smooth matching of both expansions provides a check of both. 
We clarify the discrepancy between the recent determinations of the full NNLO QCD correction to 
the semileptonic 6 — > c rate, and its earlier estimate. 



. I. INTRODUCTION 

o ■ 
o ■ 

psj ' Approximately one out of five decays of the &-quark produces a c-quark accompanied by leptons. Those semileptonic 
^ decays provide information about quark masses, the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vcb, as well 
O as properties of hadrons containing heavy quarks. In order to access that information, measurements of the decay 
probability and distributions are compared with theoretical predictions that account for radiative corrections, quark 
masses, and non-perturbative effects of strong interactions. Given that the strong coupHng constant at the mass-scale 
of the 6-quark is quite large, as {q^ = ml) ~ 0.2, and that the present uncertainty in Vcb approaches the one-percent 
level it is important to determine the second-order effects, O (a^). 

The full C (a^) correction to the decay rate was first calculated in the limit of a massless produced quark (relevant 
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I ■ for the decay b ulvi) Effects of the c-quark mass were known, until recently, only in the so called Brodsky- 
Q-i' Lepage-Mackenzie (BLM) approximation 0, [J, estimating the largest part of the second order corrections using the 
^ running of as- The remaining, non-BLM corrections, are usually smaller and much more difficult to determine. They 
/ were known only in three special points of lepton kinematics: the zero recoil, where the leptons are emitted back-to- 
back and the produced quark remains at rest 0,0,01) the maximum recoil, with the vanishing invariant mass of the 
(N ■ leptons and the intermediate recoil, where the invariant mass of the leptons equals that of the c-quark In 

^ , the latter study, the information from all three points was used to estimate the O {af) correction to the total decay 
■ rate with a massive c-quark. 

Very recently, two independent studies determined the full mass dependence of the non-BLM corrections: in the 
' calculation was performed numerically for arbitrary quark masses, and in [l^ an expansion around small rric/mi, was 
obtained analytically. The two methods are very different, with the former being more accurate at large, and the latter 
at small mj., but they agree very well in the physically interesting region of TOc = (0.25 . . .0.30) to;,. Unfortunately, 
the resulting non-BLM correction disagrees almost by a factor of two with the estimate found in 

The goal of the present paper is twofold. First, we want to check the intermediate-recoil expansion presented in [lo| . 
' Among the three kinematical points on which the estimate of the total correction was based, the intermediate- 
recoil is the only one not checked by an independent calculation. An expansion is constructed from the opposite limit 
than in [l3|: whereas there the expansion was around the zero- recoil limit, here we start from the vanishing nic, as 
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5h ■ shown in Fig. [TJ In addition, the old expansion around the zero-recoil limit is repeated and extended to higher orders. 
. . . Our second goal is to clarify the source of the disagreement between the three-point estimate and the recent explicit 
calculations. 

Fig. [T] puts the present expansion in the context of the possible kinematics of a heavy to light quark decay, Q 
q + (W* ivi). Along the diagonal, the mass of the virtual W* is equal to the light-quark mass. The arrow 
originating from the zero- recoil line corresponds to the expansion done in (and repeated in the present paper), 
while the dashed arrow coming from the zero mass point corresponds to the expansion presented here for the first 
time. Ultimately, these expansions should give a consistent value for the decay width r{b cW*). It is related to 
the differential semi-leptonic width, 

-q V{b ^ cW )\^2^w,)^^2, (1) 
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where q is the invariant mass squared of the leptons, and Fermi constant is G_f = 
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II. EXPANSION FROM ZERO MASS POINT 

Using the intermediate-recoil relation mw' — Wc, we here calculate the width as a series in p = ^ ^ 1, and as, 
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T{b^cW*) = To 



Xo + Cf—Xi + Cf(—) X2 + 0(c 



where 



9i\Vct\'mt 
647rm2(W^*)' 
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The tree-level and first-order results Xq^i are known exactly [13, llJ|, and the present approach, described below, has 
been tested with them up to 0{p^°), 



Xo = {1- p") v/1- V, 
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Xi = ^-y ( y-^-91npl +pM91np-- ) + 
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To evaluate the 0{a1) corrections, we considered the imaginary parts of 39 three-loop self-energy diagrams with 
massive propagators, such as in Fig. O and used the optical theorem to calculate the decay width. To deal with the 
two scales, mf, and rric, we used the method of asymptotic expansion , [iBl- As an example of how this asymptotic 
expansion is done, consider the left hand diagram in Fig. [51 We consider "regions" where each loop momentum is either 
hard, '-^ mt, or soft, ^ nic, and Taylor expand the propagators so that in the end we only have to deal with single 
scale diagrams as shown in Fig.[3l This method produced as many as 11 regions for a single topology. Expansions to 
the desired order 0{p^^) required algorithm [l3| for the un-factorized three loop regions, (e.g. Region 1 in Fig. [H). 

The second order results can be separated into a sum of gauge invariant parts, each with a different color factor. 



X2 = TRiNLXL + NsXs + NHXH) + CFXF + CAXA. 



(6) 



In this expression, the X^'s are the gauge invariant parts in terms of p, Nl is the number of quarks lighter than a 
c-quark, Ns and Nh serve as markers to separate the c-quark and 6-quark loop contributions. Cf = |, Ca = 3, 
and Tr — i, are the appropriate color factors in ^[/(S). The contributions from a top quark loop are not considered 
here because they are suppressed by the ratio m\,/mt and are neghgible. Terms up to C(p^°) have been calculated 
completely analytically. Here we present the formulas for terms up to to save space, Eqs. ([Il [U [9l [lOl [TT]) , while 
the numerical coefficients of all terms are given in Table [H 
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Figure 1: The kinematic region where the decay Q — > qW* is allowed. The solid arrows show known expansions while the 
dashed arrow shows the expansion presented here. The decay width is also known analytically along the whole zero recoil line. 
The dotted line corresponds to the decay width Q qlvt. The three circles along this line show known values coming from 
the different expansions. In the case considered in this paper, M = nib and m = rric. 
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Figure 2: A sample of the diagrams needed for the intermediate-recoil expansion presented here. 
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^ , ^ ■ ^ 3-2 , 


[1] = A;^ + [2] = ki + mt [3] = fe^ + 
[4] = {k'i - kif + ml [5] =ki+ml [6] = (fcs - fe)' 
[7] = (fc4 - pf [8] = (p + fcg - fc2 - fc4)' 


Region 1 

3' 2 


fci 2> rrtc, fc2 > mc, fca » mc 


[1] ^ fcf, [2] ^ [3] ^ fc| 


Region 2 

2 


fell 2> TTlc, 1^2! >• JTlc, Ifcsj ~ TTlc 


[l]^ki, [2]^kl [4]^kt+m't 
[6]^ fci = [2], [8] ^ (p-fe-fc4)' 
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|fci ~ mc, fc2 >• mc, jfcsl S> rric 


[2]^kl [3]^kl [4] ^ k'i + ml 
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' '8 ^ ^ " 


fci ~ mc, fc2 » mc, fca ~ mc 


[2] ^ ki, [4] ^ m^, [6] ^ fc2^ 
[8] fc| 



Figure 3: The asymptotic expansion of a two-scale diagram used to integrate the left-hand diagram in Fig. [2] The dashed, thin 
and thick lines correspond to massless, soft-scale massive and hard-scale massive propagators respectively. 
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For this expansion, we have used the MS definition of as normaHzed at the pole mass nib. 
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Table I: Numerical coefficients of the expansion presented here to all orders calculated. 
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III. EXPANSION FROM THE ZERO-RECOIL LINE 



An alternative way to compute at the intermediate recoil is to expand around the zero-recoil limit where rric — 
niw* — The decay width parameterization in Eqs. (|l|2p . as well as the decomposition of the second order 

correction into color parts, Eq. are still valid. For the purpose of the expansion around the zero-recoil limit, it is 
convenient to parameterize the dependence on the quark variable in terms of a new variable, (3 = 1 — 4/?^, and pull 
out its square root, thus defining new functions A^, 

X^ (p) = v//3A» (/?) , i = 0, 1, 2, A, F, L, S, H. 

The expansion around /3 = was first carried out in [1(J]. Our purpose in this section is to repeat that calculation, 
extend it to higher powers in /?, and make sure that the results match the expansion around the zero-mass point, 
p = 0, presented in Section [ill The one- loop correction in the /? expansion reads 



Ai = ^ln2-3 + /3f-ln2+iln/3-- 
8 \8 2'^ 48 

^2/^28, „ 7, ^ 13483\ /44, „ H , ^ 143263\ 

— In 2 H n /? ] + 3^ i — \n2 + — \n 3 . 

\15 15 ^ 7200 J \35 35 ^ 117600/ 

In the strong coupling constant was normalized at the geometrical mean of the quark masses, asi^/mbmc), while 
here we use as{mi,), in order to be able to match with the expansion around p — 0. Also, in [10||, the c-quark and 
6-quark loop contributions were added together and denoted Ah, while here we separate them. The 6-quark loop 
contribution is denoted by Ah and the c-quark by As. For reference, the A5 and Ah terms are given in Eqs. lfT2l[T3|) 
up to order (both normalized with Ug (mi,)), 
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where R2 is obtained from Q and has a numerical value of R2 ~ —0.72964. 
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Figure 4: A sample of the diagrams needed to calculate the expansion presented in [T3l and updated here. 
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Table II: Numerical coefficients to all orders calculated for the updated expansion from the zero-recoil limit. The values have 
been calculated using as(mi,). 

While these changes were carried out, an error was noticed in the charge renormahzation used in In that 

paper, as was normaHzed at y/m^^nTc. The error consisted in using five quark fiavors to run down to that scale, 
instead of excluding the 6-quark in the range between to;, and ^/mbm^. This error originates in {&\ . We have corrected 
for this in Eq. ^ and Table III 

To have proper matching between the expansion in and the expansion presented here, we also found that the 
old expansion needed more terms than could be obtained with the available computing resources in 1998. We have 
updated the expansion to include analytical terms up to j3^ as compared to previously. To carry out this calculation, 
we used the same methods as the authors of Instead of calculating the corrections using self-energy diagrams 
and the optical theorem, we calculated each second-order decay diagram seperately. This required the calculation of 
73 diagrams with zero, one or two loops and up to four-particle phase space integrations. Fig. |4l 

In this expansion, the loops have been integrated using the same methods described earlier. This lead to the 
calculation of 14 regions with only one having an eikonal propagator [3| ■ The numerical coefficients of all terms that 
have been calculated here are shown in Table [TTl 

Fig. [H shows how the updated expansion diflFers from the previous one and clearly displays the need for the higher 
order terms, at and below the physical value p ~ 0.3. 

In an attempt to account for the higher order terms, the authors of added a term to approximate the remainder 
of the series equal to the product of highest order term and 2(1-13) ■ This also gave an estimate of the error in their 
calculation. For a value of p = 0.3 (/3 = 0.64) they found, 

^/^2 = -4.72(14). (14) 

With the extra terms we have calculated here and the corrections to the charge renormahzation, this changes to, 

V^A2 = -4.45(1), (15) 

where we have used the same method of estimating the error. With an error of « 0.2%, we have suflflcient accuracy 
for computing the full decay width T{b ciiyg) in the next section. 

Comparing the two expansions, around p = and around p = 5, as shown in Fig. [6l one can now see that all of 
the different colour contributions and thus the full corrections agree very well. 
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IV. ESTIMATE OF THE FULL CORRECTION TO THE SEMILEPTONIC DECAY RATE 

A. Notation 



So far in this paper we have been concerned with the decay of a 6-quark into a c-quark and a virtual VF-boson, at 
the intermediate recoil where the masses of c and W* are equal. We now want to use the results we have obtained, 
together with previously obtained values at zero- and maximum- recoil to fit the corrections to the decay b civi. 
We follow the notation of fl3|, 
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In addition, we define the corrections for the integrated decay rate. 
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As we have already mentioned in the Introduction, the NNLO corrections A2 and B2 can be divided into the BLM 
and the non-BLM parts. 



A2 = Tr{NlAl+NsAs + NhAh) + CfAf + CaAa 

_ ^BLM _|_ ^nBLM 
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and similarly for the integrated corrections B. All the functions A in the Eq. iflGj) depend on two variables: the quark 



mass ratio p and the invariant mass of the leptons ^y q^. The full dependence on these variables is not yet known. The 
expansions such as described in this paper and earher studies determine A's along the sides and the bisector of the 
triangle shown in Fig. [TJ Of particular interest are their values along the vertical line corresponding to the physical 
value of p ~ 0.3, describing the differential decay rate dr{b civf)/dq^ . Ref. [l3| used the three known points along 
that fine to fit a polynomial and integrate Eq. (U) over q^ , thus providing an estimate of the second order non-BLM 
corrections to the full semi-leptonic decay width. 



BLM 



0.9(3). (19) 

This numerical value is quoted from [ll|, where the author fixed another mistake in (lol| related to using the maximum- 
recoil result. Recently, however, two calculations of the full decay width T{b cive) [ll|, gave 



BLM 



1.73(4). 



(20) 




Figure 5: Second order corrections expanded from the zero recoil line. The dashed line shows the expansion up to 0{f3'^) while 
the solid line shows the expansion up to 0{(3^). For the purpose of comparing with this plot is made assuming Qs(ym(,mc) 
is used in the NLO correction. 
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This section is devoted to clarifying the discrepancy between these results. 

For comparison purposes we use as{^Jmb^mZ), Nf = 4, and p = 0.3 as ii i lioll . where Nf refers to the number of 
light quarks used to calculate the BLM contributions. The authors of [ill, Es] use a different set of parameters, to 
which we will return in section BVDI 



B. Effect of corrected coupling constant normalization 

A large part of the discrepancy between Eqs. I|19|20p is due to the incorrect charge renormalization, as discussed 
in Section IIIIl We have corrected this and recalculated the non-BLM contributions using the same fitting method 
described in [lO|. Analogously to Eq. (8) in [10], we introduce a new function of the lepton invariant mass at fixed 




0.1 0.2 0.3 0.4 0.5 ^ 

Quark loop contributions 




Total 0(0^) correction. 

Figure 6: The matching between the different colour contributions and total X2 contribution. The thick line corresponds to 
the expansions presented here and the thin line corresponds to the updated expansion from the zero recoil line. 
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quark-mass ratio p (we use p = 0.3). It is denoted ^{q ) and defined as 



2, A2(g2)-AfLM(g2) 



^Born(g^) 

The three available values of S,{q^) at = 0, m^, and q^^^^ — {mi, — mc)^ are 

e(0) = 1.26, ami) = 1.27, ^('zLx) = 0.19. (22) 

These values have been obtained using results from and 01, with the &-quark charge renormalization terms 

from corrected. Fitting these values to a function defined by, 

C(<z') = aig-* + azg' + 03, (23) 

we integrate over q^ to find a value for the non-BLM corrections. The values quoted in Eq. I|22p are normalized to 
the Born rate, ^Bom, so the integral needed is analogous to Eq. (9) in [lo| . 

^nBLM ^ /o"^^^rfg'^Born(g')e(g') ^24) 
* /o'"'"'d<72^Born(g2) 

This integration, with the input from Eq. gives B"^^^^ = 1.4(2). This agrees with Eq. |[20]) much better than the 
value given in Eq. lfT9|) . The error is estimated by performing the analogous fit of the BLM corrections and comparing 
the result to the exact value 



C. Effect of extending the expansion to higher orders in /9 

In section HV BI we have merely corrected the renormalization in the old results. In addition, using the results of the 
two expansions in the present paper, we can obtain a more accurate input along the intermediate recoil line. Instead 
of the value 1.27 in Eq. (|22|) . this gives ^(m^) = 1.33. This is related to the shift induced by the higher-order terms, 
illustrated in Fig.[5l There we see that the full correction is less negative than previously assumed, thus the difference 
with the BLM correction is more positive (larger). Since the zero- and the intermediate-recoil points are close to each 
other, even a small shift of the value at one of them may be amplified in the integral of the fitted function. 

After the integration in Eq. l(24|) . this change leads to the new value B^f^^ — 1.5(2) which now agrees with 
Eq. I|20p . The error estimated by comparing with the BLM case is about 12 per cent. By correcting the mistake in 
renormalization and including more terms in the expansion from zero recoil, we have brought the disagreement from 
about a factor of two down to w 10 per cent, within the quoted error margins. 



D. A better fitting method 



Further improvement is possible using a better method of fitting the polynomial. In Eq. (|2T|) . we normalized the 
points to the tree level result Asoin- We find that, if this normalization is not done, i.e. instead we use, 

aq') = A,iq')-A^^^\q% (25) 

the polynomial fit gives a much better estimate of the exact result. With this adjustment of the fitting procedure, we 
end up with a final non-BLM estimate of iJj^^^ = 1.76(4), a significant improvement. Without knowing the shape 
of the dT{b clvi)/dq^ curve, we cannot say whether this is a numerical coincidence. We have also performed this 
fitting for the 0{as) corrections and BLM approximation. Both estimates give results that are within « 3 per cent 
of the exact known result and are more accurate than using an analog of Eq. l(2T|) for the fit. 

In the more recent papers Jill, [l3| the authors use a different set of parameters for their calculations, namely 
as(mb), Nf — 3, and p — 0.25. For easy comparison, we have also calculated the non-BLM corrections with this 
set of parameters. Using Eq. (|25| for the fitting procedure and the expansion about p = presented here, we find 
^nBLM ^ 3.37(15)^ as compared with the result of B"BLM ^ 3 40(7) from 
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V. SUMMARY 



To summarize: we have corrected an error in the strong coupling constant normalization in the previous 
intermediate-recoil expansion. We have extended that expansion to several higher orders in the parameter /3, de- 
scribing the difference between nic and mb/2. We have confirmed the correctness of that expansion by constructing 
a new one, also along the intermediate-recoil diagonal but around its other end, corresponding to mc/mi, 0. This 
analysis allowed us to re-evaluate the fit of the dr{b civi)/dq^ curve based on the three kinematical points, and 
remove the disagreement between the correction to the total rate T{b cii'e) obtained from such a fit, and that 
obtained from the direct four-loop calculations (ill , [l^ . With this result, the full massive calculation of the O (a^) 
corrections to the semileptonic 6-quark decay rate is confirmed. 
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Appendix A: CONTRIBUTIONS FROM b cccW* 



The expansions used to obtain the maximum recoil point for our polynomial fit were calculated in Q and The 
two expansions agree very well except for the region with TOc < ^ ■ This discrepancy can be attributed to the omission 
of the amplitude of 6 ^ cccW* from the expansion in 0. For completeness, we have calculated this contribution for 
both maximum recoil and intermediate recoil, as this expansion was also not included in These expansions are 
calculated as threshold expansions in terms of S given by, 

mc = ^(l-.5), (Al) 

where the 3(4) indicates the factor used when calculatin g th e maximum recoil, 3, or intermediate recoil, 4, expansion. 
The methods used for both expansions are discussed in [ia|- This calculation relies on the ability to reduce the four 
particle phase space integrals needed, into a product of two particle phase spaces. 

For the maximum recoil case, the expansion has been calculated up to (5^^, with the first four terms given here, 

r[b^cccW )m^.=o = ^ (A2) 

83 , 7 ,2 55 ^3 753 ,4 

H S H 6^ H S'^ H 6^ + . 

56 64 896 896 

For intermediate recoil, the expansion has been calculated up to S^\/S, 

535 , 137045 175277863 

5 5^ H 5^ + 

108 85536 13343616 
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